Gorkov equations for a pseudo-gapped high temperature superconductor 
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A theory of superconductivity based on the two-body Cooperon propagator is presented. This 
theory takes the form of a modified Gorkov equation for the Green's function and allows one to 
model the effect of local superconducting correlations and long range phase fluctuations on the 
spectral properties of high temperature superconductors, both above and below Tc- A model is 
proposed for the Cooperon propagator, which provides a simple physical picture of the pseudo- 
gap phenomenon, as well as new insights into the doping dependence of the spectral properties. 
Numerical calculations of the density of states and spectral functions based on this model are also 
presented, and compared with the experimental STM and ARPES data. It is found, in particular, 
that the sharpness of the peaks in the density of states is related to the strength and the range 
of the superconducting correlations and that the apparent pseudo-gap in STM and ARPES can be 
different, although the underlying model is the same. 

PACS numbers; 74.20.-z, 74.25.-q 



I. INTRODUCTION 



The anomalous properties of high temperature super- 
conductors (HTS) have been the object of many investi- 
gations over the last ten years B Among these anomalous 
properties, the pseudo-gap phenomenon (a substantial 
decrease of the one particle density of states near the 
Fermi energy in the normal state below a certain tem- 
perature T*) has been studied by several experimental 
techniques. In tunneling spectroscopy experiments, in 
particular, the pseudo-sfp manifests itself as a gap in 
the excitation spectrumJj On the theoretical front, many 
competing models have been proposed.tltBl One of the 
popular interpretations of the pseudo-gap is that super- 
conductivity forms locally at T* , but the phases of distant 
superconducting "droplets" remain incoherent until 
superconducting transition temperature Tc is reached 
This view is supported by an increasing evidence that 
the pseudo-gap phenomenon is intimately connected to 
the underlying superconducting phase, mainly because 
the d-wave symmetry of the pseudo-gap is the-same as 
that of the gap in the superconducting phase.ll3llil The 
phase fluctuation model of the pseudo-gap state is dif- 
ferent from the usual theory of superconducting fluctu- 
ations, because the latter is only valid Jifiar Tc, and in- 
volves both size and phase fluctuations .Ej 

A theory of the pseudo-gap above Tc has also conse- 
quences below Tc- In particular, it is connected to the 
character of the excitations responsible for destroying su- 
perconductivity. Different mechanisms (thermal phase 
fluctuations, quantum phase fluctuations, nodal quasi- 
particles) may all contribute to the properties in the su- 
perconducting state, and these contributions may be of 
varying importance if pae considers low temperatures or 
temperatures near Tct3 One may add that the clue to 
a theory of high temperature superconductivity will go 
through the explanation of detailed properties, like the 



absence otfluasiparticles above Tc and their appearance 
below Tcj^l or the_anomalous properties of the density of 



states in vorticesii3 

Within the BCS theory, the inclusion of phase fluctua- 
tions in a droplet model must be done in two steps. First, 
one introduces a local BCS gap, with a phase, and then 
one must somehow analyze the phase fluctuations in a 
second analytical siep distinct from the first. The theory 
of Franz and Millialj gives an example of this type of ap- 
proach. Starting from the form of the Green's function 
in a uniform superflow, they extend it semi-classically to 
non-uniform situations assuming slow spatial variations 
of the superfluid velocity. This Green's function is then 
averaged over a Gaussian distribution of velocity fluctua- 
tions, which relates the result to a correlation function of 
the velocities.. A similar approach has been used by Kwon 
and Dorsey,c3 who treat the coupling to the fluctuating 
phase using a self-consistent pexturbation theorj^As em- 
phasized by Geshkenbein et alLj and Randeria,El strong 
pairing correlations should however be incorporated at 
the basis of any model of the pseudo-gap regime. 

A systematic theory of the effect of phase fluctuations 
on the density of states (and other properties) in HTS, 
above and below Tc, must start by putting the phase- 
phase correlation function at the core of the theory for 
superconductivity, at the same level as the size of the 
local gap. This program implies that one avoids devel- 
oping the theory of superconductivity by defining a gap 
function and an anomalous propagator. This is equiva- 
lent to writing the BCS theory in a particle number con- 
serving scheme {i.e. in the canonical ensemble, avoiding 
the definition of an anomalous amplitude between states 
with different particle numbers). This theory has actu- 
ally bepi written down forty years ago by Kadanoff and 
MartinlH3 (KM), and rediscovered by others, in particular 
for the discussion of Josephson arrays.E3 The KM theory 
is based on the two-body Cooperon propagator, and de- 
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scribes quite naturally the effect of phase fluctuations. 
This theory has already been appliect-to the HTS in a 
nice series of papers by Levin et al.f3 but with a dif- 
ferent interpretation (not related to phase fluctuations), 
a different focus, and a different formalism than in our 
work. In this paper we rewrite the basic KM equations 
for the case of a lattice Hamiltonian, we show again how 
the standard BCS theory is recovered with a straightfor- 
ward approximation for the two-body Cooperon correla- 
tion function, and we then derive the basic equations for 
a pseudo-gap state, in particular the equivalent Gorkov 
equations which have to be used in the calculation of 
vortex states. 



II. KADANOFF-MARTIN EQUATIONS IN A 
LATTICE MODEL 

We consider the lattice Hamiltonian 

ijiT ij 

In this, c|g. creates an electron at site i and &J = c-j^cj| 
creates a Cooper pair at site i. We assume that t and V 
are symmetric and real. The usual Gorkov equations can 
then be written as a single equation: 

[g°M] '^GM = 1 + EK) g{co„), (2) 

where Gijir) = -(TV{Cj|(T)c]^(0)}} and 

= - E ^^„,s.,i?;,v;.,a° (-^„). (3) 

r-ir-2 

[S^l^n)] ^^."'^ = iWnSij — Uj is the free Green's function, 
Wn are the odd Matsubara frequencies, and = (bi). 
The notation [G'~'] G implies matrix multiplication in 
the {i, j} space. These equations are supplemented by 
the self-consistent equation for Bi (or A^) and the self- 
consistent equation relating the number of particles N to 
the chemical potential. The equation for Bi is 

with.F*(r) = -(T,{ct^(T)c]^(0)}). 

The Kadanoff-Martin correlation function description 
of superconductivity consists simply (after a long and 
thorough discussion of higher order correlation functions) 
in replacing S in Eq. (j|) by 

^tj{uJn) = -^^^V,^^L^^^^{uJn + UJm)V^^jGji{uJm) 

(5) 



where Lr2ri{T) = (^V {^^ (''')^r2 C-*)}) Cooperon 
propagator. The self-consistent equation for B is re- 
placed by the equations 

E E + i-^n)V^^L,r, i^rn) (6a) 

for T>Tc, and 

~^ E E ^^r2 + ^n)Gt2 Vi (^™) (^b) 

for T < Tc, where flm are the even Matsubara frequen- 
cies. The fact that the inhomogeneous term in the lad- 
der equation for L has to be dropped when calculating 
the order parameter in the condensed state is not much 
commented upon in the original paper of KM, but it is 
related to the range of Lr2rii^m,), which is finite above 
Tc and infinite below Tc (see below). In fact, the quantity 
J2rir2 ^r2ri{^m) IS of the Order N"^ below Tc (where N 
is the number of sites) whereas the corresponding sum of 
the inhomogeneous term in Eq. ( |6a| ) is only of order N. 
The usual BCS theory is recovered by setting 

Lr2ri{T) = B,._^B*^, Lr2rii^m) = PB.j.^B*^5n^fi 

in Eq. (H). The self-consistent equation Eq. ( |6b| ) for L 
then goes over into the self-consistent equation Eq. (^) 
for B, and clearly Eq. (||) goes into Eq. (|). In the BCS 
framework, the Thouless criterion for Tc becomes the gap 
equation for T < Tc- 

The KM description of superconductivity, which is en- 
tirely based on the properties of the function L, is thus 
seen to be a natural starting point if one wants to in- 
troduce explicitly local order and phase fluctuations in 
the physical description of high temperature supercon- 
ductors. This is done in the next section. 



III. PHENOMENOLOGICAL DESCRIPTION OF 
A PSEUDO-GAPPED SUPERCONDUCTOR 

Our fundamental assumption is that Eq. (||) is gener- 
ally valid, in the sense that it expresses in general the 
single particle Green's function in terms of the Cooperon 
propagator, regardless of the model or the approxima- 
tions involved in calculating this propagator. Our pur- 
pose in this paper is to explore the experimental conse- 
quences of a simple heuristic form for L, which is the 
translation of the physical picture presented in the Intro- 
duction. For T > Tc, we write: 

Lr2rAr) - |BO|2i^(rl2/^^o) + \B^? F{r2 - n, r) (7a) 
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with ri2 = \r2 — Ti] and R{x) some cutoff function which 
vanishes rapidly for x > \. This equation expresses the 
fact that there are strong superconducting correlations 
at the scale go, represented by a finite value of S", going 
to zero gradually at a temperature T* . The strength of 
the superconducting correlations between "droplets" is 
represented by an amplitude and some function F 
describing essentially the correlations in an XY model 
above the Kosterlitz-Thouless (KT) transition. Both 
and F will be temperature dependent. As T approaches 
Tj, from above (we identify with the KT transition 
temperature) , the correlation length of F diverges and F 
approaches 1. For T < Tc, we write: 



(7b) 



where only the amplitude B^ carries a phase. The 
assumption here is that short range correlations have 
a strong incoherent part, even in the superconducting 
state. When introduced into the equation for the self- 
energy, Eq. (7b) for L means that the self-energy in the 
superconducting state is the sum of a coherent and an in- 
coherent part; this appears to be the case in sotue recent 
calculations based on a fermion-boson model. Inspec- 
tion of Eqs. (7b) and (^ shows that if, in the ladder 
approximation for a homogeneous system, it turns out 
that Lr2ri is the sum of a constant term and a term of 
finite range, then the constant term will obey Eq. (|6b|), 
whereas the finite range term will obey Eq. (|6a|). 
Eqs. (0) translate into the following equation for Q: 



(8) 



where 



+ J2 V.r, R{ri2/go) K,,g%{-CJn) (9) 



and 



E E - ri, u;„ + io^) 

xV,^^g°,iiUrn), T>T, (10a) 

-Y,V,,^BlXX.,5U-^n), T<T,. (10b) 



Inspection of Eq. (10b), together with Eq. (g), shows 
that below Tc, a pseudo-gapped superconductor obeys 
the following modified Gorkov equations: 

{ [g'M] ''gM\ , +E VuBjT* (i^r.) = 6,, (lla) 

{[c;°(-c.„)]''^*(c^n)}..-E^^^^r^.,K) = 0. (lib) 

Quantum Monte Carlo (QMC) calculations of the pair- 
ing correlations were recently reported for the attractive 



Hubbard model at zero temperature.c3 Although these 
results are restricted to short distances (~ 6-10 lattice 
sites) we tentatively connect our model to the QMC cal- 
culations with the following arguments. For the sys- 
tem sizes considered in the QMC calculations (typically 
14 X 14 sites) the correlation fmiction for the largest dis- 
tance in the system has not converged to its asymptotic 
value. We attribute the slow decrease of the correlations 
at intermediate distances (see inset of Fig. 4 in Ref. |2^) 
to a large value of qq with respect to the system size. The 
results of Ref. ^ also show that the strength of the pair- 
ing correlations at intermediate distances increases, and 
differs increasingly from the BCS result, as the Hubbard 
interaction U/t increases. Closer inspection of the data 
in Fig. 3 of Ref. |26| indicates that the ratio of the BCS to 
the QMC correlations at intermediate distances is also an 
increasing function of U/t. The simplest BCS approxi- 
mation to Eq. ( |7b| ) is to replace the cutoff function R 
by 1 , describing correlations which are independent of ri 
and r2 (in a homogeneous system, the second term of 
Eq. (7b) is a constant |S^P). With this approximation, 
one can account for the above trends by assuming that 
both B'^ and the ratio B^ /B^ increase as U/t increases. 
Finally, we shall include in our numerical calculations the 
"onsite" correlations found in Ref. ^ for distances within 
two lattice spacings, by adding a term 



Klr, = {\Br + me 



-2ri2/a 



(12) 



to the model Eqs. (|^), where a is the lattice parameter. 
We find, however, that this correction has a negligible 
impact on the spectral functions, and could equally be 
dropped without changing the results presented below. 



IV. NUMERICAL RESULTS 

We now use the general equations derived in Section 0, 
together with the model Eqs. (Q) and (^2|), to calculate 
the temperature dependence of the density of states and 
spectral functions in a homogeneous system. The calcu- 
lations are compared with the STM and ARPES experi- 
mental results for BSCCO. In order to reduce the num- 
ber of adjustable parameters, we take for the correlation 
function F{r, r) — exp[— r/^(T)], which describes time- 
independent phase-phase correlations above Tc. The cor- 
relation length ^ is equal to Qiae-xjp\b/ yL{T — Tc)/ J], 
with a w 0.21, b « 1.73, and J « rc/0-89.E3 The length- 
scale Qi is the lattice parameter of the effective 2D ~ XY 
model describing phase fluctuations. We found that the 
main features in the spectral properties above Tc are 
rather insensitive to the details of the correlation func- 
tion. We have explicitly checked this point by comparing 
different functions F^ in particular functions which give 
a better description of the correlations in the XY model. 
The cutoff function R is modeled as exp(— r/gio)- 

The four parameters po, Qi, and B^ are chosen to 
achieve good agreement with the experimental results. 
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S(fc, UJn) 



V\q)L(q)dq 



g — fc 



(13) 



FIG. 1. Model two-body Cooperon correlation function 
at temperatures below and above Tc. The increase of the 
correlations for r < a is due to L"'^ given in Eq. (p^. Below 
Tc, L{r) converges to the finite asymptotic value \B^^ at 
distances of the order go if > and of the order a if 
B° = 0. Above Tc, the range of L{r) is finite. If B° > 0, 
this range is given by max{go, ^(T)) while if _B" vanishes it is 
given by ^(T). 



We use the value qq — 50a; if the first term in the right- 
hand side of Eq. ( [7b| ) is the dominant one (S*^ > B^), 
we found that a relatively large value of go is needed in 
order to obtain well developed coherence peaks in the 
zero temperature density of states. In addition, we see 
that, according to the previous discussion, go must be 
large with respect to 14a. The parameter gi controls 
the temperature evolution of the spectral functions above 
Tc, and takes the value gi = 5a. The larger gi, the wider 
the temperature region above Tc in which finite range 
phase coherence contributes to the pseudo-gap. The am- 
plitude is adjusted to fix the gap energy to ~ 40 meV. 
Finally, the ratio B'^ / B^ is varied in order to control the 
relative importance of short range superconducting corre- 
lations and long-range phase fluctuations. The behavior 
of the resulting Cooperon propagator L{r) is illustrated 
in Fig. for temperatures below and above Tc and for 
different values of \B°/B^\. 

For a translationally invariant system and our model 
Cooperon propagator, Eq. (||) can be recast as 



where V{q) = Vo + 2Vi{cosqxa + cosqya), L{q) is the 
Fourier transform of L{r), and is the free dispersion. 
Here, Vo and Vi are the onsite and nearest-neighbor po- 
tentials, respectively, and we neglect next-nearest neigh- 
bor interactions; we assume Vi = Vo/4 in all of our cal- 
culations. For the dispersion, we use a tight-binding ex- 
pression which fits the BSCCOj-Eermi surface and corre- 
sponds to a bandwidth of 2 eVo The self-energy at real 
frequencies is evaluated by making the analytic continua- 
tion iujn — *■ =0^-1- zO^ in Eq. (O) and discretizing the 
BZ integralo The spectral function is then calculated ac- 
cording to A{k, Lj) = — ilm{[a;+ — — S(fe, ct;+)] 
and the density of states is N{uj) oc Jg^ ^(^i ^) dk. It 
is easy to check from Eq. (^3|) that, if L{q) > — a 
condition obeyed by our model — then the Green's func- 
tion is analytic in the upper half of the complex plane, 
the spectral function A{k, u) is positive, and the Green's 
fimction goes to zero as lo^^ for \uj\ — > oo. 



A. Scanning tunneling spectroscopy 

Neglecting possible anisotropics of the tunneling ma- 
trix element as well as fc^-dispersion effects, we calculate 
the tunneling conductance as the convolution of the den- 
sity of states with the derivative of the Fermi function. 
The result is shown in Fig. ^ for various temperatures. 
In order to focus on the effect of local superconductivity 
and phase fluctuations, we have kept the model parame- 
ters independent of temperature: the whole temperature 
dependence of the curves, in Fig. |, relates to the varia- 
tion of the correlation length ^ and Fermi function with 
T. A better flt to the experimental data could be ob- 
tained, in principle, by allowing the amplitudes B'^ and 
B^ to vary with temperature. This would not, however, 
change the qualitative conclusions we wish to draw. In 
Fig. |(a), B° is larger than B^ while in Fig. |(b) B^ 
is larger than B^ . In the next section, we argue that 
these two typical cases correspond to underdoped (UD) 
and overdoped (OD) situations, respectively. The spec- 
tra shown in Fig. ^ reproduce some of the characteristic 
features observed experimentally in BSCCO samples.cl 
Both UD and OD curves evolve smoothly across Tc into 
a pseudo-gapped spectrum, the peak-to-peak distance re- 
maining approximately temperature independent. More- 
over, the coherence peaks and the gap structure disap- 
pear more rapidly in the OD case as the temperature 
is raised, which is also consistent with the experimen- 
tal flndings. The model, however, is not able to account 
for a number of experimental observations, such as the 
asymmetry in the temperature dependence of the pos- 
itive and negative-bias conductance peaks, or the dip 
structure recorded at ~ 2A below Tc- We also note that 
the model Eqs. (Q) has s-wave symmetry. The calculated 
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FIG. 2. Tunneling conductance as a function of tem- 
perature. The model parameters represent (a) under- 
doped (VoB° = 15 meV, B° / = 2) and (b) overdoped 
(VoB" = 7 meV, B°/B^ = 0.5) situations. The critical tem- 
perature is Tc = 80 K (bold line). The dashed lines show 
the T = 300 K spectra corresponding to =0 and (a) 
VoB° = 7.5 meV, (b) VbS" = 3.5 meV. The curves have been 
shifted for clarity. 



spectra are therefore not expected to agree in details with 
experiment at low energies. 

According to our model, the local superconducting cor- 
relations responsible for the high temperature pseudo- 
gap also have implications below T^. In the underdoped 
case, the local (incoherent) superconducting correlations 
broaden the zero temperature density of states. The re- 
sulting conductance spectra have small coherence peaks 
and a rounded line-shape around the fermi energy. In the 
overdoped case, in contrast, the T = curve looks more 
like a (s-wave) BCS spectrum. 
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FIG. 3. Calculated zero-bias conductance as a function of 
temperature for underdoped (UD, black symbols), and over- 
doped (OD, white symbols) systems. The parameters are the 
same as in Fig. ^. The dashed lines show the conductance 
obtained by thermally broadening the T = Q spectra. The 
dotted line is obtained by letting B^ go to zero in the OD 
situation. The vertical axis starts at zero. 



As the temperature increases from zero to Tc, the den- 
sity of states remains unchanged in both UD and OD 
cases, and the temperature dependence of the conduc- 
tance spectra relates solely to the Fermi function. This 
behavior persists above Tc in the UD case owing to the 
dominant role of (which is T-independent in our cal- 
culations). In the OD situation, on the contrary, the 
gap fills in rapidly above Tc as the contribution of -B^ 
disappears due to increasing phase fluctuations; at ele- 
vated temperatures, only a weak pseudo-gap due to 5° 
remains. Fig. || also illustrates the effect of the tem- 
perature dependence of the amplitudes and B^ . At 
room temperature, B^ is expected to vanish and i?'^ is 
expected to be smaller than at low temperature. Tak- 
ing fii = and S°(300 K) = S"(0 K)/2, one obtains 
the dashed spectra in Fig. 0, which no longer exhibit a 
sizable pseudo-gap structure. 

The difference between the temperature evolutions of 
the UD and OD spectra is best seen in Fig. |, where we 
plot the calculated zero-bias conductance. Below Tc, the 
zero-bias conductance is larger in the UD case due to 
strong local correlations. Above Tc, the conductance in- 
creases sharply in the OD case, corresponding to the fill- 
ing of the gap. In either UD and OD cases, the zero-bias 
conductance above Tc is larger than the value expected 
by thermally broadening the T = spectra (see Fig. H). 

From a general point of view, one can confirm from our 
calculations that the sharpness of the peaks in the den- 
sity of states (and correspondingly the size of the zero- 
bias conductance) is related to the strength and range 
of the superconducting correlations. The larger the ra- 
tio B^ /B'^ and/or the longer the range go, the sharper 
the peaks (the smaller the zero-bias conductance). As 
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FIG. 4. Calculated ARPES intensity near (vr, 0) as a 
function of temperature for underdoped (UD) and overdoped 
(OD) systems. The parameters are the same as in Fig. ^ 
The curves have been shifted for clarity. Inset: representa- 
tion of the Brillouin zone showing the Fermi surface used in 
the calculations and the Fermi crossing near (tt, 0). 



an example, we show in Fig. ^ the conductance obtained 
by letting the coherence term go to zero in the OD 
situation. Comparison of the curves with and without 
shows that the phase coherence has the effect to de- 
press the density of states at the Fermi energy — there- 
fore raising the coherence peaks — below Tc and in some 
temperature range above T^, where the correlation length 
C(T) is large. 



B. Angle-resolved photoemission 

Experimentally, it is found that the temperature de- 
pendence of the energy dispersion curves measured by 
ARPES near (tt, 0) depends on doping. In overdoped 
samples, the leading-edge midpoint energy moves toward 
the Fermi energy — suggesting that the gap closes — 
as T increases above Tc- The temperature variations of 
the mid po i nt energy are usually smaller in underdoped 
samplesMN 

Apart from a matrix element, the ARPES intensity 
is just the product of the spectral and Fermi functions. 
This quantity, calculated at the Fermi crossing near the 
M point, is shown in Fig. || as a function of tempera- 
ture. A clear difference between the temperature evolu- 
tion of the spectral peak in the UD and OD cases can 
be seen. Consistently with experiment, the peak shifts 
toward the Fermi energy in the OD case as T increases. 
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FIG. 5. Temperature dependence of the quasiparticle en- 
ergy (maximum of the curves in Fig. ^ at the Fermi crossing 
near (tt, 0) for underdoped (UD, black circles), and overdoped 
(OD, white circles) systems. The squares show the energy of 
the quasiparticle at a fe point just below the Fermi surface 
along the AI~Y line. 



In the UD case, the peak position is to first approxima- 
tion independent of temperature. Below Tc, the curves 
are almost identical to the spectrum at Tc (because the 
temperature T < Tc « 7 meV is small with respect to 
the peak energy 35 meV) and are not shown. One can 
see that the quasiparticle peak is much sharper at Tc in 
the overdoped as compared to the underdoped system. 
This has also been seen experimentallyEll and can eas- 
ily be understood in our model. The destruction of long 
range order by phase fluctuations clearly affects qualita- 
tively the spectral functions in the OD case where the 
transition across Tc is accompanied by a decrease of the 
quasiparticle lifetime and increase of the intensity at the 
Fermi energy. 

The position of the main quasiparticle peak in Fig. ^ is 
reported in Fig. ^ as a function of temperature. The tem- 
perature dependence of the gap was studied in Refs |l^ 
and |3l| by fitting the experimental ARPES curves to a 
three parameters Green's function. For overdoped sam- 
ples, the gap was found to decrease with increasing tem- 
perature (Ref. |3^), in a way very similar to what we 
obtain in the OD case, although the decrease was found 
to begin already below Tc- Note that a small finite gap 
persists at all temperatures in our calculations, since no 
temperature dependence of B'^ and B^ was taken into ac- 
count. In a real situation, B'^ and B^ would both vanish 
at some temperature above Tc- In the underdoped sam- 
ples, the gap was found to be temperature independent 
within error bars (Ref. or slightly decreasing above 
Tc (Ref. The trend in Fig. |is similar. The sHght 
decrease of the gap above Tc in the UD case results from 
the suppression of the phase correlations as the temper- 
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ature is raised. 

The spectral line-shapes in Fig. || are considerably 
sharper than what is usually measured by ARPES, es- 
pecially at elevated temperatures. The estimated ex- 
perimental resolution of ^ 10 meV cannot alone explain 
this difference. Similar conclusions have been reached in 
Ref. It was shown there that the inverse quasiparticle 
lifetime implied by fitting the experimental spectra are 
an order of magnitude larger in ARPES with respect to 
STM. Inhomogeneities-in the sample properties could ex- 
plain this discrepancy,L3 since a much larger region of the 
sample surface is probed by ARPES compared to STM. 

The temperature dependence of the OD quasiparticle 
peak in Fig. ^ contrasts with the apparent temperature 
independence of the gap width in Fig. ^. The coherence 
peaks in the density of states are due to quasiparticle 
states with momenta just nearby fcp- Therefore, one may 
expect that the energies of all these quasiparticles evolve 
in the same way as the temperature increases. In this 
case, the coherence peaks would rigidly follow this tem- 
perature evolution and both the STM and ARPES gaps 
would close in the OD situation. We have found, how- 
ever, that in our model the energies of the quasiparticles 
at and nearby kp have different temperature dependen- 
cies in the OD case. This is illustrated in Fig. ^, where 
we plot the energy of a quasiparticle with a momentum 
k just below the Fermi surface along the M-Y line. For 
T < Tc, this particular fe point contributes to the co- 
herence peaks in the density of states, since the corre- 
sponding energy is within 2 meV of the energy at kp. 
At 200 K, instead, the two energies differ by ~ 5 meV 
in the OD case, which is approximately half the width 
of the zero temperature coherence peaks. This explains 
why the STM gap fills in instead of closing although the 
ARPES gap at fcp closes. Thus, our results show that the 
apparent "visual" pseudo-gap may be different in STM 
and ARPES data, even if each measurement is in agree- 
ment with the same underlying theory. 



V. CONCLUSION 

Many workers in the field share the belief that the 
pseudo-gap phase in HTS is a kind of mixed state, where 
strong short range superconducting correlations coexist 
with long range phase disorder. This fact should reflect in 
the properties of the Cooperon propagator, which should 
show "partial" superconductivity even above Tc- In this 
paper, we have shown that it is possible to describe the 
properties of pscudo-gapped superconductors by writing 
the superconductivity theory in general in terms of this 
Cooperon propagator, and that reasonable phenomeno- 
logical assumptions about the form of this propagator 
lead to good agreement with experimental data. We have 
thus a theoretical framework which is valid both above 
and below without special treatment of the pseudo- 
gapped phase. Our main assumption is that the relevant 



difference between overdoped and underdoped HTS is in 
the relative magnitude of the short and long range parts 
of the Cooperon propagator, described by the parame- 
ters 5° and B^, respectively. In the underdoped HTS 
we assume that the ratio B'^/B^ is larger than in the 
overdoped HTS. We tentatively claim that is related 
to the single-particle energy gap Ap measured by single- 
particle spectroscopy, while B^ is related to the coherence 
gap Ac measured in Andreev reflectipn or Josephson ex- 
periments. As shown by DeutscherJ^ Ap and Ac differ in 
the HTS: the ratio Ap/Ac is close to one in the overdoped 
region and increases as the doping is reduced. 

In this paper, we do not attempt to calculate the 
Cooperon propagator using one or the other theoretical 
method. We first want to derive some empirical con- 
straints on the function L from direct comparison with 
experiments. Our approach is also limited, at this stage, 
to s-wave gap symmetry. We are currently working on 
an extension of these calculations for d-wave symmetry 
and on the calculation of the density of states in vortices. 
Also, comparisons of our model with other detailed spec- 
troscopic data below Tc (vortices and Josephson effect 
in particular) will show whether the approach presented 
here is a fruitful one. 



ACKNOWLEDGMENTS 

We wish to thank S. E. Barnes, H. Beck, 0. Fischer, 
M. Franz, B. W. Hoogenboom, and J.-M. Triscone for 
very useful discussions. 



* Bernard . Giovannini@physics . unige .ch 
^ For a recent review see, e.g., T. Timusk and B. Statt, Rep. 

Prog. Phys. 62, 61 (1999). 
^ Ch. Renner, B. Revaz, J.-Y. Genoud, K. Kadowaki, and 0. 

Fischer, Phys. Rev. Lett. 80, 149 (1998). 
^E. Roddick and D. Stroud, Phys. Rev. Lett. 74, 1430 

(1995). 

" V. J. Emery and S. A. Kivelson, Nature 374, 434 (1995). 
^ P. W. Anderson, The Theory of Superconductivity in the 

High-Tc Cuprates (Princeton University Press, Princeton, 

1997). 

^ M. Randeria, J.-M. Duan, and L.-Y. Shieh, Phys. Rev. 
Lett. 62, 981 (1989). 

J. R. Schrieffer and A. P. Kampf, J. Phys. Chem. Solids 

56, 1673 (1995). 
* A. V. Chubukov, D. Pines, and B. P. Stojkovic, J. Phys.: 

Condens. Matter 8, 10017 (1996). 
^ P. A. Lee and X.-G. Wen, Phys. Rev. Lett. 78, 4111 (1997). 

J. Ranninger, J. M. Robin, and M. Eschrig, Phys. Rev. 

Lett. 74, 4027 (1995), and references therein. 

V. B. Geshkenbein, L. B. loffe, and A. I. Larkin, Phys. Rev. 

B 55, 3173 (1997). 



7 



A. Perali, C. Castella ni, C. Pi Castro, M Grilli, E. Piegari, 
and A. A. Varlamov, cond-mat/9912363. 

" A. G. Loeser, Z.-X. Shen, D. S. Dessau, D. S. Marshall, 
C. H. Park, P. Fournier, and A. Kapitulnik, Science 273, 
325 (1996); J. M. Harris, Z.-X. Shen, P. J. White, D. S. 
Marshall, M. C. Schabel, J. N. Eckstein, and I. Bozovic, 
Phys. Rev. B 54, R15665 (1996). 

" H. Ding, T. Yokoya, J. C. Campuzano, T. Takahashi, M. 
Randeria, M. R. Norman, T. Mochiku, K. Kadowaki, and 
J. Giapintzakis, Nature 382, 51 (1996). 
For some recent efforts to separate the effect of phase fluc- 
tuations from size fluctuations, see P. Curty and H. Beck, 
Phys. Rev. Lett. 85, 796 (2000). 

A. Paramekanti and M. Randeria, ;ond-mat / 0001 10g| . 
A. Kaminski, J. Mesot, H. Fretwell, J. C. Campuzano, M. 
R. Norman, M. Randeria, H. Ding, T. Sato, T. Takahashi, 
T. Mochiku, K. Kadowaki, and H. Hoechst, Phys. Rev. 
Lett. 84, 1788 (2000). 

Ch. Renner, B. Revaz, K. Kadowaki, L Maggio-Aprile, and 
0. Fischer, Phys. Rev. Lett. 80, 3606 (1998). 

1^ M. Franz and A. J. Millis, Phys. Rev. B 58, 14572 (1998). 

^°H.-J. Kwon and A. T. Dorsey, Phys. Rev. B 59, 6438 
(1999). 

21 M. Randeria, cond.mat/9710223. 

22 L. P. Kadanoff and P. C. Martin, Phys. Rev. 124, 670 
(1961). 

2"^ L. Weiss and B. Giovannini, Helvetica Physica Acta 55, 
468 (1982). 

2^* B. Janko, J. Maly, and K. Levin, Phys. Rev. B 56, R11407 
(1997); Q. Chen, I. Kosztin, B. Janko, and K. Levin, Phys. 
Rev. Lett. 81, 4708 (199 8); Q. Chen, I. Kosztin, and K. 



Levin, 



:ond-mat/9908362; J. Maly, B. Janko, and K. Levin, 



Phys. Rev. B 59, 1 354 (1999); L Kosz tin, Q. Chen, Y.-J. 
Kao, and K. Levin, cond-mat / 9906180 ; K. Levin, Q. Chen, 



;ond-mat/0003133. 



and L Kosztin, 

25 S. E. Barnes, Int. J. of Mod. Phys. B 13, 3478 (1999). 

26 ^ 



M. Guerrero, G. Ortiz, and J. E. Gubernatis, cond- 



mat/9912114 . 

^' R. Gupta and C. F. Baillie, Phys. Rev. B 45, 2883 (1992). 

2* J. Schmalian, S. Grabowski, and K. H. Bennemann, Phys. 
Rev. B 56, R509 (1997). 

2^ The self-energy Eq. ( |l3| ) is most easily calculated using a 
Fast Fourier Transform algorithm. To achieve a good pre- 
cision, we set 0^ = 1 meV and use a 1024 x 1024 fe-point 
mesh in the Brillouin zone for calculating the density of 
states, and a 2048 x 2048 mesh when calculating individual 
spectral functions. 

^'^ M. R. Norman, H. Ding, M. Randeria, J. C. Campuzano, 
T. Yokoya, T. Takeuchi, T. Takahashi, T. Mochiku, K. 
Kadowaki, P. Guptasarma, and D. G. Hinks, Nature 392, 
157 (1998); J. Phys. Chem. Solids 59, 1888 (1998). 

^1 M. R. Norman, M. Randeria, H. Ding, and J. C. Cam- 
puzano, Phys. Rev. B 57, R11093 (1998). 

^2 M. Franz, private communication. 
G. Deutscher, Nature 397, 410 (1999). 



8 



